SU(4) Skyrmions and Activation Energy Anomaly 
in Bilayer Quantum Hall Systems 
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The bilayer QH system has four energy levels in the lowest Landau level, corresponding to the 
layer and spin degrees of freedom. We investigate the system in the regime where all four levels 
are nearly degenerate and equally active. The underlying group structure is SU(4). At v = 1 the 
QH state is a charge-transferable state between the two layers and the SU(4) isospin coherence 
develops spontaneously. Quasiparticles are isospin textures to be identified with SU(4) skyrmions. 
The skyrmion energy consists of the Coulomb energy, the Zeeman energy and the pseudo-Zeeman 
energy. The Coulomb energy consists of the self-energy, the capacitance energy and the exchange 
energy. At the balanced point only pseudospins are excited unless the tunneling gap is too large. 
Then, the SU(4) skyrmion evolves continuously from the pseudospin-skyrmion limit into the spin- 
skyrmion limit as the system is transformed from the balanced point to the monolayer point by 
controlling the bias voltage. Our theoretical result explains quite well the experimental data due 
to Murphy et al. and Sawada et al. on the activation energy anomaly induced by applying parallel 
magnetic field. 
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I. INTRODUCTION 

Exchange Coulomb interactions play key roles in vari- 
ous strongly correlated electron systems. They are essen- 
tial also in quantum Hall (QH) systemsU0|, where spins 
are polarized spontaneously as in ferromagnets. Much 
more interesting phenomena occur in bilayer QH systems 
due to exchange interactions. For instance, an anomalous 
tunneling current has been observed between the two 
layers at the zero bias voltage. It may well be a man- 
ifestation of the Josephson-like phenomena predicted a 
decade agoQ. They occur due to quantum coherence 
developed spontaneously across the layers 01 

Electrons in a plane perform cyclotron motion un- 
der strong magnetic field B± and create Landau levels. 
When one Landau level is filled up the system becomes 
incompressible, leading to the integer QH effect. In QH 
systems the electron position is specified solely by the 
guiding center X = (X, Y) subject to the noncommuta- 
tive relation, [X, Y] = —i£%, with the magnetic length 
£b = y/h/eB±. The QH system provides us with a sim- 
ple realization of noncommutative geometry. It follows 
from this relation that each electron occupies an area 
2ttI\ labelled by the Landau-site index. Electrons be- 
have as if they were on lattice sites, among which ex- 
change interactions operate. 

In QH ferromagnets charged excitations are spin tex- 
tures called skyrmi ons Ml ■ Skyrmions are identified 
experimentally^, 0- bv measuring the number of 
flipped spins per one quasiparticle. Indeed, by tilting 
samples, the activation energy increases by the Zeeman 
effect, which is roughly proportional to the number of 
flipped spins. On the contrary, an entirely o ppo site be- 
havior has been observed by Murphy et al.[ll| in the 
bilayer QH system at the filling factor v — 1, where the 
activation energy decreases rapidly by tilting samples. It 
is called the activation energy anomaly because of this 



unexpected behavior. Note that we expect an increase 
since the v = 1 bilayer QH system is a QH ferromagnet 
with spins spontaneously polarized. This anomalous de- 
crease has been argued to occur due to the loss of the 
exchange energy of a pseudospin texture based on the 
bimeron modelffj], where the spin degree of freedom is 
frozen. Here, the layer degree of freedom is referred to 
as the pseudospin: It is said to be up (down) when an 
electron is in the front (back) layer. 

We investigate physics taking place in the lowest Lan- 
dau level (LLL). Since each Landau site can accommo- 
date four electrons with the spin and layer degrees of 
freedom, the underlying group structure is enlarged to 
SU(4). Let us call it the isospin SU(4) in contrast to the 
spin SU(2) and the pseudospin SU(2). We study charged 
excitations in the v = 1 bilayer QH system. A natural 
candidate is the SU(4) isospin texture to be identified 
with the SU(4) skyrmion 0,0]. A specific feature of the 
SU(4) isospin texture is that it is reduced to the spin and 
pseudospin textures in certain limits. 

At v = 1, electrons are transferable between the two 
layers continuously without breaking the QH effect [l3|. 
Namely, the bilayer QH system can be continuously 
brought into the monolayer QH system by changing the 
density imbalance. It implies that a pseudospin texture 
at the balanced point must be continuously transformed 
into a spin texture at the monolayer point. It is natural 
that a pseudospin texture evolves into an isospin texture 
and then regresses to a spin texture in this process. The 
corresponding continuous transformation of the activa- 
tion energy has already been observed experimentally by 
Sawada et al.[l3l ll3.|l*5ll. In this paper we explain these 
experimental data |lll UM Il5| based on the excitation of 
SU(4) skyrmions. 

In Section II we derive the Landau-site Hamiltonian 
governing bilayer QH systems, by extending the alge- 
braic method employed previously to investigate multi- 
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component monolayer QH systems [l(fl. In Section III we 
rewrite the Landau-site Hamiltonian into the exchange- 
interaction form. The Coulomb potentials associated 
with the direct and exchange interactions are obtained 
analytically. 

In Section IV the ground state is explored in the 
regime where the SU(4) isospin coherence develops spon- 
taneously. It is shown that the capacitance energy con- 
sists of two terms arising from the direct interaction, 
which is the standard one made of two planes, and from 
the exchange interaction, by way of which the capaci- 
tance energy becomes considerably smaller than the stan- 
dard one for a small layer separation. 

In Section V the excitation energy of an electron-hole 
pair is calculated exactly. A hole (electron) is the small 
size limit of a skyrmion (antiskyrmion), which is real- 
ized when the Zecman gap Az and the tunneling gap 
Asas are very large. The exchange energy is very differ- 
ent whether the pair excites the spin or the pseudospin. 
Due to this difference the pseudospin excitation occurs 
at the balanced point even if the tunneling gap is quite 
large. Our result explains why pseudospin flips were ob- 
served in the experiment due to Terasawa et al.|l5j| in a 
sample having a very large tunneling gap (Asas — 33K). 
In Section VI we study skyrmions in a microscopic the- 
ory, extending the approach|l7j previously known for the 
monolayer system. 

In Section VII the effective Hamiltonian is derived from 
the Landau-site Hamiltonian by making the derivative 
expansion. We introduce composite bosons together with 
the CP 3 field^l to describe coherent excitations. In Sec- 
tion VIII we carry out an analysis of SU(4) isospin tex- 
tures identified with topological excitations called CP 3 
skyrmions or equivalently SU(4) skyrmions. A spin tex- 
ture (spin-skyrmion) and a pseudospin texture (ppin- 
skyrmion) are special limits of an SU(4) skyrmion. We 
estimate the excitation energy of one SU(4) skyrmion. 
We show that, if a ppin-skyrmion is excited at the bal- 
anced point, it evolves continuously into a spin-skyrmion 
at the monolayer point via an SU(4) skyrmion as the 
density imbalance is increased between the two layers. 

In Section IX we investigate the activation energy 
anomaly by applying the parallel magnetic field between 
the two layers. We calculate explicitly the loss of ex- 
change energy of one SU(4) skyrmion, which is shown to 
be proportional to the capacitance term. Our theoretical 
result explains quite well both experimental data due to 
Murphy et al.[lj and Sawada et al. [TllT^ . 

Section X is devoted to discussions. 



II. LANDAU-SITE HAMILTONIAN 

Electrons perform cyclotron motion under perpendic- 
ular magnetic field B± . The number of flux quanta pass- 
ing through the system is = B±S/$v, where S is the 
area and <!>d = 2irh/e is the flux quantum. There are 
N$, Landau sites per one Landau level, each of which is 



associated with one flux quantum and occupies an area 
S/N& = 2-k1 2 b . In the bilayer system an electron has two 
types of indices, the layer index (f, b) and the spin index 
(Ti !)• One Landau site may accommodate four electrons. 
The filling factor is v = with N the total number 

of electrons. 

The kinetic Hamiltonian is 

Hk = ^M J * f W( Dx ~ iD v)(P* + ^y)*( x )' 

(2.1) 

where \& stands for the four-component electron field, 
and Dk = —ihdk + &Ak is the covariant momentum. The 
kinetic Hamiltonian, which is invariant under the global 
SU(4) transformation, creates Landau levels. Assuming a 
large Landau level separation we focus on physics taking 
place in the LLL. Then we may neglect the kinetic energy 
since it is common to all states. 

We introduce 4x4 matrices r^ pm and rP pin 
generating the spin and pseudospin (ppin) SU(2) 
groups, as summarized in Appendix A; in particular, 
T s / in =diag(l, -1,1,-1) and T^' m =diag(l, 1,-1, -1). 
To avoid confusions we use S a for the spin SU(2) field, 
P a for the pseudospin SU(2) field (a — x, y, z) and I a for 
the isospin SU(4) field (A = 1, 2, • • • 15). 

The total Hamiltonian consists of the Coulomb term, 
the Zeeman term, the tunneling term and the bias term. 
The role of the bias term is to transfer electrons from 
one layer to the other by applying the bias voltage Vbias 
between the two layers. 

The Coulomb interaction is decomposed into the 
SU(4)-invariant and SU(4)-noninvariant terms, 

H C=IJ d 2 xd 2 y V+(x - y)p(x)p(y), (2.2) 

=2 fd 2 xd 2 yV-(x-y)P z (x.)P z (y), (2.3) 

where p(x) = #t(x)#(x), 2P a (x) = #t( x ) T PPi*,$r ( x ) and 

V ± M = —( , 1 A (2-4) 

1 ' Sirs [\x\ ^/\ X \ 2 +d 2 J ^ ' 

with the layer separation d. 
The Zeeman term is 

Hz = -A z 5,(x), (2.5) 

where A z is the Zeeman gap and 2S a = *V| pin *. The 
tunneling and bias terms are combined into the pseudo- 
Zeeman term, 

H PZ - -Asas^(x) - eV bias P,(x), (2.6) 

where Asas is the tunneling gap. 
The total Hamiltonian is 

H = H+ + Hc+H z + H PZ . (2.7) 
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We investigate the regime where the SU(4)-invariant where 
Coulomb term Hq dominates all other interactions. Note 
that the SU(4)-noninvariant Coulomb term vanishes, 
Hq — > 0, in the limit d — > 0. 

We expand the electron field operator by a complete 
set of one-body wave functions <^j(x) = (x|z) in the LLL, 



JY,!, 



(2.8) 



where c a {i) is the annihilation operator at the Landau 
site \i) with a =ff, fj, bt and b|, 

{c a (i),cj 3 (j)} = 5 i:j 5 a p, 

{c a (i), C/3 (i)} = {ct(i),4(i)}=0. (2.9) 

It is impossible to choose an orthonormal complete set 
of one-body wave functions <^(x) in the LLL|lj|- Hence 
-0 Q (x) does not satisfy the standard canonical anticom- 
mutation relation, as implies that an electron cannot be 
localized to a point within the LLL. Various interactions 
are projected to the LLL 0,0, EJH(| by expanding the 
electron field as in J2|S}. 

The projected density is given by p(x) = '!>''' (x)'I'(x) 
with the use of the field operator (|2.8|l . Its Fourier trans- 
formation reads 161 



= * -4q 2 /4 



m\e ^In)^^, 



(2.10) 



where X = [X, Y) is the guiding center obeying the non- 
commutativity [X, Y] — —ii 2 B . Similar formulas hold for 
spin density operators and so on. Substituting them into 
(|2.2() and (|2.3() we find the projected Coulomb terms to 
be 

H£ = E y «U E 4(m)4(z) C/3 (j> a (n), (2.11) 

mnij a (3 

h c = E v^na E i^ p ' m u[rr n u 

mnij 



x ci(m)cUi)c s (j)cf3{n) 



(2.12) 



with 



V 



III II I J 



= 4^ J d 2 pV±(p)(m\e^\n)me-^\j) (2.13) 



and 



1±e -|p|d 



?P / 2 



(2.14) 



8tt£ |p| 
We also find 

flz=-A z £)S , ,(m,ro), (2.15) 

JTi 

m 

(2.16) 



p(m,n) = ^2cl(m)c a (n), 

a 

Sa(m,n) J- £ 4(m)[rf n ]^ C/5 (n), (2.17) 



a/3 



and a similar formula for P a (m,n). These formulas are 
derived precisely in the same way as for the multicom- 
ponent monolayer system with the replacement of the 
potential V mn ij by V^ njj -: See Section V in Ref.[l^|. 

For a later convenience we represent the projected 
Coulomb Hamiltonians H2.11fl and (|2.12l) as 



mnij 

H D =4 E V- nij P z (m,n)P z (i,j). 



(2.18) 
(2.19) 



III II i J 



(2.20) 



In the momentum space the total Coulomb Hamiltonian 
reads 

H U =n J d 2 P V+(p)p(-p) P (p), 

+ 47T Jd 2 pV B (p)P z (-p)P Z (p), 

where we have defined 

1 oo 

ftp) =^E( m i e " ipX i n )^ m '")' 

mn 
1 oo 

^(p) =^EHe~ ipX w^K™)- 



(2.21) 



We call (|2.20|1 the direct-interaction form of the Coulomb 
Hamiltonian He- 
Note the relation 



p(q)= e -^ 2 / 4 p(q). 



(2.22) 



between the two types of densities (|2.10|) and 12.21fl . 
Though /o(p) presents a useful tool to work with noncom- 
mutative geometry 16], it turns out that p(q) describes 
the physical density [see l|5.19|l ] . 

The density operator p(p) satisfies the Woo algebra[l8l 

MM, 



[?(p),?(q)] = -p(p + q)sin 



iPAq 



(2.23) 



For the SU(4) isospin field we may define 

I A {m,n) = ^4HM«(J^W, (2-24) 

a/3 

where A a is the generating matrices [See Appendix A], 
and 



4 



All of the density operators /5(p) and Ia(p) form an al- 
gebra, which we have called 16j Woo(4) since it is the 
SU(4) extension of Woo. It gives the fundamental alge- 
braic structure of the noncommutative system made of 4- 
component electrons. In general there arises the Woo(N) 
algebra in the N-component QH ferromagnet. 

We define classical densities p cl (p), I A l (p), S^'(p) and 
PM by 



p cl ( P )He|/Kp)|6), i c 2 = (6\i a (p)\&) 



(2.26) 



and so on, where |6) represents a skyrmion state. In the 
coordinate space the relation 



1 



^ /l 1 (x)*^(x) = ^? C1 (x)* 

A=l 



N 
2%t% 



p w 



(2.27) 

holds among the classical densities associated with the 
generators of the Woo(N) algebra, where *k stands for 
the Moyal star product. We demonstrate this formula in 
Appendix [B] Remark the invariance of (|2~2T|) under 



P cl (r) 



N 
2ir£l 



(2.28) 



This represents the electron-hole symmetry. 



III. EXCHANGE INTERACTIONS 

In classical theory the Coulomb energy is simply given 
by (|2.2I) with the use of the classical density /9 cl (x), but 
this is not the case in quantum theory. The exchange 
interaction emerges as an important interaction from the 
exchange integral over wave functions. We present a rig- 
orous treatment of the direct and exchange energies valid 
in QH systems. 

For this purpose we rewrite the microscopic Coulomb 
Hamiltonians 1)2. and l|2.12ll into entirely different 
forms. Based on an algebraic relation, 

1 N2 ~ X 1 

SafjSra = 77 ^ [Aa] ot [^a] a(3 + J^S aT d a p, (3.1) 
A 

which holds for SU(N) with the generating matrices, 
these Coulomb Hamiltonians are equivalent to 



15 



- 2 E F "U [ E WaV, n) 



■mm j 



A=l 



+ 2jfP(m,j)p(i,n)], (3.2) 



and 



15 



H x = - 2 E V ™nij [ E UlA(m,j)I A (i 



mnij 



A=l 



+ ^y^( to ' 3)p(h n)] , (3.3) 



where N = 4, Ij\(m,n) is defined by (|2.24() and 

£1,2,3,8,13,14,15 = +1) S4,5,6,7,9,io,ii,i2 = — bi tnc mo- 
mentum space they read 

„ 15 

H+=ir / c1 2 pV+(p)[YIa(-p)Ia(p) 

J A=l 

+ 2^?(-p)p(p)], (3.4) 

. 15 

H- =4n / d 2 P V x (p) [ ]T UIa(-p)Ta(p) 

J A=l 



+ 2]V?(-P)?(P)], (3-5) 



where 



y±(p) = ^L Jd 2 ke- u zP Al <V±(k) (3.6) 
and N = 4. 

We change the SU(4) basis from \ A to rf in , rPP in and 
T^ pin T^ pm by way of the formula (|A6|) in Appendix, which 
transforms variable Ja(to,7i) to a set of variables, 

Sa(m,n) = i^ct(m)(rf n ).r C rW, 

err 

P a (m,n) = ^4(m)(rPP in ) ffr c r (n), 

err 

R ab (m, n) = ±Yl cUrn)(rT n rf m Ucr(n). (3.7) 

The sum of l|3.2() and (|3.3|l is expressed as 

Hx = - E V mnij[ S a(m,j)S a (i,n) + P a (m,j)P a (i,n) 

mnij 

+ Rab(rn,j)Rab(i,n)] 
~ 2 E V m~nij [S a (m,j)S a (i,n) + P z (m, j)P z (i, n) 

mnij 

+ Raz{m,j)R az {i,n)} 
~ \ E v mnijP{mJ)p(i,n), (3.8) 

mnij 

where 

V ■ ■ = V - * - -4- V~ V d = V" >r — V~ f3 Q"! 

v mmj v mnij ' v mnij 5 v mnij v mnij v mnij' V "^/ 

In the momentum space we find 

H X = - I / dW(P)[£a(-p)£a(p) + P (-p)P„(p) 

+ J R ab (-p) J R a6 (p)] 

-7T / d 2 pF x (p)[S Q (- P )S a (p) + P 2 (-p)P z (p) 



+ i? az (-p)i? a2 (p)] 

rf 2 p^x(p)?(-p)p(p), 



(3.10) 
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where 

Vx (P) - V+ (p) + v- (p) , vi (p) = V+ (p) - v- (p) . 

(3.11) 

In EHJ and EHJ the summation (Ea= xyz >E b = xyz ) 
over the repeated indices a and 6 is understood. We 
call H3.10|) the exchange-interaction form of the Coulomb 
Hamiltonian He- 

We have demonstrated that the Landau-site Hamilto- 
nian Hq possesses two entirely different forms, the direct- 
interaction form Ho given by (|2.2U|) and the exchange- 
interaction form Hx given by QiS.lUp . They are equiva- 
lent, Hd = Hx, as the microscopic Hamiltonian. In this 
paper we are interested in the excitation energy (6\H\6) 
of a skyrmion state |S). In a previous paper [l6| we have 
presented a heuristic argument showing that 



ie\H a \e) = m 



m 



(3.12) 



Here H$ and HQ are the Hamiltonians of the direct- 
interaction form l|2.2U[l and of the exchange-interaction 
form (|3.1U|) . where the density operators p(p), Ia(p), 
■ ■ ■ are replaced by the classical ones p cl (p), Ia(p)> 
We present a proof of this formula for the skyrmion state 
in Appendix iBl Thus, the energy of a charge excitation 
consists of two well-separated pieces, the direct energy 
and the exchange energy. 



IV. GROUND STATE 

We first determine the ground state of the bilayer sys- 
tem in this section. For simplicity we start with the spin- 
frozen SU(2) bilayer system. We are concerned about 
the regime dominated by the SU(2)-invariant Coulomb 
interaction Hq. Hence we define the ground state as 
an eigenstate of ff<t given by l|2.11|l . and treat all other 
interactions as small perturbations. 

The unperturbed ground-state energy is given by 



with 



1 



jnnj I 



(4.1) 



1 / <P*v£(k) 
1 [i ±e d 2 /24 erfc ^/V2^) 



— rv ± . 

jn 



A 



where 



A 



2 Att££ b 



(4.2) 



(4.3) 



We take A^, as the Coulomb energy unit in this paper. 
The unperturbed system is equivalent to the monolayer 
SU(2) QH ferromagnet, where all pseudospins are spon- 
taneously polarized into one arbitrary direction in the 



SU(2) pseudospin space. The ground state may be ex- 
pressed as 



-'- Wo/ V^4(») 10), (4-4) 



by introducing two constant parameters <7o and t?o- We 
call (To the imbalance parameter since it represents the 
density imbalance between the two layers. The ground 
state is a coherent state due to an infinite degeneracy 
with respect to Co and 

We then study the effect due to the SU(2)-noninvariant 
interaction. First, due to the capacitance effect (oc 
V-mnij) au pseudospins are polarized into one arbitrary 
direction within the pseudospin xy plane, implying that 
the electron density is balanced between the two layers 
(do = 0). The symmetry SU(2) is broken into U(l). It is 
thus said that the bilayer QH system is an easy plane fer- 
romagnet. Next, we apply the bias voltage to generate an 
imbalanced density state (co 7^ 0), where the symmetry 
is still U(l). Finally, the tunneling interaction breaks the 
symmetry completely by fixing $0 = in the state (|4.4|l . 
Consequently, the ground state is the bonding state |B) 
parametrized by the imbalance parameter ao, which is 

dHU) with #0 = Q- 

The generalization to the SU(4) bilayer system is triv- 
ial by adding the spin component. The ground state 
turns out to be the up-spin bonding state |B|) due to the 
Zeeman effect. It is convenient to represent the up-spin 
bonding state |B|) and the up-spin antibonding state 
|Af) as 

|BT)=n B T(")l°)' |AT)=rK(")l°>' ( 4 - 5 ) 



where 



B T (n) \ = J_ ( VTTTrTt 



Co 



y l - go 



c fT (n) 



(4.6) 

The down-spin bonding state |B|) and the down-spin an- 
tibonding state |A|) are similarly defined. 

We evaluate the ground state energy E g = (B1 \H\B^\) 
by representing the Hamiltonians in terms of these oper- 
ator. The ground-state energy per one site reads 



EJN* 



r^cap^Q 



-AgAsV 1 -^- ^biasCTo- Az, (4.7) 



where 



with 



4(e D - £ x) 



£n = 



> v~ •■ 

nnjj 

jn 



1 d 



2tt It 



A 



(4.8) 



(4.9) 
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and e x given by l|4.2[) . The imbalance parameter <jq is 
determined to minimize the ground-state energy, 

AsAS^O 



eVbia 



(4.10) 



MS i\ function of the bias voltage. 

The second term in the ground-state energy (|4.7J) is 
the capacitance energy, 



pcap 

g 



4 2C- (4 ' n) 

where Q = ^ep ao is the charge imbalance and C is the 
capacitance per unit area. We rewrite (|4 . 1 Ul) as 

Q 
c 

with 

Co 



(4.12) 



(4.13) 



When the tunneling interaction is absent (Asas — ¥ 0), 
(|4.12l) is reduced to 

Q = CV hias - (4.14) 

Eqs. 1)4. 11|) and (|4.14|l are the well-known formulas for 
the condenser, where the bias voltage Vbi as is balanced 
with the electric potential Q/C due to the charge differ- 
ence Q between the two layers. In this case, a charge 
transfer between the two layers makes no work, since 
there exists no potential difference between the two lay- 
ers. 

When Asas 7^ 0, the cancellation is imperfect as in 
(|4.12l) . Thus, when a charge is moved from one layer to 
the other, it is necessary to supply an energy against the 
potential difference Vj unc = Vbi as — Q/C. Consequently, 
the pseudo-Zeeman energy is given by 

Hpz = - [Asas^(x) + eVj unc P 2 (x)] 



CO 



*SAS, 



H 



PZ 



E 



P x (m,m) 



CO 



=P z (m, m) 



(4.15) 



iSAS- 



(4.16) 

As Asas - ► 0, the pseudo-Zeeman energy vanishes even 
in imbalanced configuration (ctq =/= 0), and the total 
ground-state energy consists solely of the capacitance en- 
ergy H4.11f> . In order to analyze charge excitations it is 
necessary to use (|4.16H as the pseudo-Zeeman term in the 
total Hamiltonian (|2.7(l . 



V. ELECTRON AND HOLE EXCITATIONS 

As we see in the succeeding section, a skyrmion and 
an antiskyrmion are reduced to a hole and an electron 
in their small size limit. We analyze electron and hole 
excitations in this section to derive some exact results. 



r-level 
p-level 
s-leve! 
g-level 



A | 



FIG. 1: The LLL contains four energy levels corresponding 
to the two layers and the two spin states. (We call them the 
g-level, s-level, p-level and r-level.) At v = 1 the ground state 
is the up-spin bonding state. An electron may be moved to 
any one of the other three levels to form an electron-hole pair 
excitation. 



A. Electron-Hole Pair Excitation 

One electron may be excited from the up-spin bond- 
ing state into the down-spin bonding state, the up-spin 
antibonding state or the down-spin antibonding state 
[FIG^. These electron-hole states are 

IV.) = ;hf) = B{Cj)B T (*)|BT), (5.1) 

|V P ) = |ef;hf)=A{(i)S T (fe)|BT), (5.2) 

|V r ) = |ef i ;hf)=4(i)fl T (*)|Bt). (5.3) 

We assume that an electron and a hole are separated far 
enough so that interactions between them are neglected. 
The energy matrix looks as 



<V r |ff|Vs> «v|tf|Vv> o 
o o (Vp|flWp> 



(5.4) 



where 



(^ItflVO =£ g + 2e+ + 2^ex+A z , 
(^ r |if|V r ) =E g + 2e+ - 2a 2 e x - a 2 e cap + A 2 



^SAS 



A 



SAS 



(i> p \H\^ p ) =E g + 2e+ - 2e x - a 2 e c&p -r 

V 1 - CT 6 

mH\i> s ) =2a 0V /l-a2 £D . (5.5) 

We have (^> s |if|^ p ) = (^} T \H\^ V ) = 0, because the 
Hamiltonian does not flip the spin (it does not involve 
S x and S y ). 

The matrix l|5.4|l is diagonal at the balanced point. 
The minimum eigenvalue is the spin-excitation energy 



(i> s \H^ s ) = E s + 2e+ + A z , 



(5.6) 



or the ppin excitation energy 



(Vpl^p) =P g + 2e+-2e x + A S AS. (5.7) 
The pseudospin excitation occurs when 

2e x >A SA s-A z . (5.8) 
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/'b=1 1 .5 nm. d=1 1 .Snm 
d SAS =30K 




C B =1 1.5nm, d=3.1nm 




o.a 0.4 o.e o.e ? 
Density Difference $jjg 



o.2 o.4 o.e o.a i 

Density Difference Oq 



given by 

l S +> =c\(0)l[c\(n)\0), |6_) = ]Jc\(n)\0), 

(5.11) 

where we have placed an electron or a hole at the 
momentum-zero state. Their classical densities are 

p c i(m ) n)=^(6 ± |4(m)c (7 (n)|6 ± ) 

a 

=S mn ± SmoSnQ- (5.12) 

In the coordinate space, based on formula (|2.21() it reads 



FIG. 2: The energy of an electron-hole pair state is depicted 
as a function of the imbalance parameter ao based on the 
formulas ( I5.5B and 15.101 with typical sample parameters as 
indicated. Three curves corresponds to the state \ip B ), IV'sr) 
and IV'p)- We have normalized the activation energy of the 
pseudospin excitation to zero at ao = 0. 



We remark that 2e x ~ 47K in a typical sample with 
d ~ 231 nm. 

The matrix (|5.4|l is diagonal also at the monolayer 
point, where the spin excitation always occurs with the 
excitation energy 

(iP s \H\i> s )=E g + 2(e+ + ex)+A z , (5.9) 

since the energies of the other two modes diverge. 

The two excitation states |^ s ) and \tp T ) mix to make 
a new state, \ip ST ), to lower the excitation energy except 
for ao = and uq = 1. After diagonalization it is 



(1>JH hO = E s + 2e+ - ^e cap 



is 



SAS 



- - V(Asas - 4a§6 D )2 + 16^(1 - al){^f 



(5.10) 



with Asas = Asas/ v 1 — Oq< Thus, when the condition 
(|5.8|) is satisfied, an electron is excited to the \if) p ) state 
at the balanced point, and it makes a sudden transition 
to the \ip ST ) state, and finally transforms smoothly into 
the jV's) state, as depicted in FIG|2 The sudden transi- 
tion would be smoothed out in a skyrmion-antiskyrmion 
excitation involving several electrons simultaneously. 



B. Physical Densities 

It is interesting to investigate a single electron excita- 
tion and a single hole excitation separately. For simplic- 
ity we analyze the SU(2) QH ferromagnet. One-electron 
excited state |6+) and one-hole excited state |S_) are 



p c i(x)=p (l±2e- r2 



(5.13) 



where p is the electron density in the ground state. The 
hole density becomes negative at the origin, 



limpl(x) 

x^0 



-Pa- 



(5.14) 



There is nothing wrong with this mathematically since 
the electron cannot be localized within the LLL. Never- 
theless, we cannot accept this as a physical quantity. 

We recall that the wave function of one electron with 
the angular-momentum zero in the LLL is 



1 



which leads to the density 

<5p cl (x) = |^ (x)| 2 = (0o e-'- 2 /^ 



(5.15) 



(5.16) 



at v = 2ttP b p = 1. When we remove one electron from 
or add one electron to the filled up-spin level, the density 
becomes 



pg(x)= Po ±r5p cl (x), 

and the hole density satisfies 



limp! 1 (x)=0. 

x^0 



(5.17) 



(5.18) 



This behavior is what we expect for the physical density 
at the origin. 

It is easy to see that the Fourier transformations of 
(|5.13|) and l|5.17|l are related as 



Ph(q) 



(5.19) 



This is precisely the relation l|2.22|l between the two types 
of the densities. We conclude that the projected density 
represents a physical quantity. 
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C. Direct and Exchange Energies 



A. SU(2) Skyrmions 



For simplicity we work still in the SU(2) QH ferromag- 
net. The electron- hole pair excitation energy isp^| 



6E pair = (i> s \H c \^ s 



E g = 2e x , 



where 



ex = Vbi 



A 



(5.20) 



(5.21) 



as is obtained by taking d — > in 1(4.2(1 and 1)5.5(1 . It fol- 
lows from (|2. 13|) that Voooo is the direct integral and Vjooj 
(j 0) is the exchange integral. Hence, the Coulomb en- 
ergy consists of the direct energy SH^ = Voooo an d the 
exchange energy 5H$ = J2jjto V j00j- 

We can express them in familiar forms. First, the 
classical energy associated with the density modulation 
(ETT2T) reads 



6Hg = ^ 



V mnij Sp c \m,n)Sp c \i,j) = Vo o (5.22) 



with Sp cl (m, n) — 5 m oS n o, which is transformed into 

8H$ = l -J d 2 xd 2 y V(x - y)5p c \x)5p c \y) (5.23) 

with <5/9 cl (x) being given by 1(5.16(1 . This is the direct 
energy of the excitation. We next remark that the spin 
reads 

S c x \m, n) =S c y \m, n) = 0, 

Sf(m, n) =- (5 mn - S m0 S n o) (5.24) 

both for the electron or hole excitation. Substituting 
them into the exchange energy 1)3.2(1 . or 

Hi = - 2 £ V mmj [ J2 Sf(m,j)S°% n) 

ci—xyz 

+ ]pi(m,j)pi(i,n)], (5.25) 



mmj a—xyz 



we reproduce SH^ = YljM^jOOj- I n the momentum 
space we have a more familiar expression, 

H*=n I d 2 P V x (p)[ S*(-P)S?(P) 

a—xyz 

+ ^(-p)p d (p)]. (5.26) 

Thus, the excitation energy of an electron or a hole 
presents us the simplest example of the decomposition 
formula 1(3.12(1 into the direct and exchange energies. 



VI. SKYRMIONS IN MICROSCOPIC THEORY 

We study skyrmions in a microscopic theory of the 
SU(2) QH ferromagnet, and then extend the scheme to 
the SU(4) bilayer QH ferromagnet. We set lg = 1 
throughout in this section, where p = l/2n. 



A skyrmion and an antiskyrmion are topological soli- 
tons in the nonlinear sigma model, and characterized by 
the nonlinear sigma field (normalized spin field), 



ax 



S v — 



• 2 + a 2 



z ~ 2 r 2 + 



when it carries topological charge ±1. Actual skyrmions 
are much more complicated in the QH ferromagnet be- 
cause the spin rotation modulates the electron density 
due to the Woo(2) algebra. Nevertheless, a gross feature 
remains as it is[l|. 

Let 1 6) be the skyrmion state. Its spin field is given 
by [see (E23J] 



n)(G\S a (m,n)\6). (6.2) 



Due to the formula 
(m|e ikX |m + n) 



y/(m + n)\ V V2 



ky -\- ik x 



-k 2 /4rn 



(6.3) 



we need to have (&\S Xt y(rn,n)\&) oc 5 m , n ±i and 
(&\S z (m,n)\&) oc S m , n to be consistent with 1(6. l|l . Such 
a state is uniquely constructed as|l7j 



i 6sky )= n^wio) 



(6.4) 



with 



£ f (n) = u{n)c\{n) + v{n)c\(n + 1). (6.5) 



They satisfy the standard canonical commutation 
relations, 



[£(m),£(n)] = 6 mn , E(m),£(n)] = 0, 



(6.6) 



provided u 2 (n) + v 2 (n) — 1. Note the state |G s k y ) de- 
scribes one hole state in the momentum-zero site when 
we set u(n) = and v(n) = 1 for all n. 

Similarly, the antiskymion state is given bv[l7) 



|6 asky ) = nC t (n)4(0)|0) 



with 



C>) 



-u(n)c|(n + 1) + u(n)cj (n). 



(6.7) 



(6.8) 



The state |S as k y ) describes one electron excited state 
in the momentum-zero site when we set u(n) = and 
v(n) — 1 for all n. In what follow we only discuss the 
skyrmion case explicitly. 
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First we evaluate p cl (m,n) = (&\p(m,n)\&) and oth- 
ers. The only nonvanishing components are [see H-B8I) in 
Appendix [E] 

/^(m, to) = u 2 (m) + v 2 (m — 1), 

Sf(m, to) = — — [u 2 (to) — v 2 (m — 1)] , 

S^(to, m + 1) = 5^(to + 1, m) = -iS%\m, m + 1) 

= iS^Om + 1, to) = —u(m)v(m), (6.9) 

where we have set v(— 1) = 0. They are converted into 
the momentum space based on the formula 12. 211 . 

? cl (k) =^ e - fc2/4 E KM + « 2 (™ - !)] £ ™ ( fc2 /2) , 

m=0 

S?(k) = - -le- fc2 / 4 ]T [« 2 (to) - ^ 2 (to - 1)] L m (fc 2 /2) 

m=0 

2*00 =^e~*'* £ !^)^) ( fc 2 /2) , 
2v27T ^— 1 vm + 1 

v m— v 

Sf (k) = ^ e -^/4 V <^MiW (fc 2 /2) (6.10) 

v m— v 

after some calculation with the use of (16.31) . 

The electron number dcffcrcncc between the skyrmion 
state and the ground state is 



where we have set x = (r cos 6, r sin 8). The physical 
densities reads 



1 



p cl (x) =—e 



-r 2 /2 



S£(x) 



S?(x) 



EK(-) + - 2 (-- 1 )](y) 

m=0 

2tt ^ 



u 2 (m) — v 2 {m — 1) /r 2 
2~ 



1 -r 2 /2 u(m)v(m) /r 2N 



2vrV2 



n\\Jn + 1 



.14) 



<5Q = E [p cl (™,m-) - 1] 

m=0 

oo 

=u 2 (0) - 1 + E [« 2 (to) - m 2 (to - 1)] = -1 



(6.11) 

by pairwise cancellations, where we have used v 2 (m) — 
1 — u 2 (m), to > 0, and v(— 1) = 0. Hence the skyrmion 
excitation removes one electron from the ground state. 
However, because the manipulation is too subtle, we 
give a concrete calculation taking an explicit example 
in (HITHJ- 

As we have remarked before, it is necessary to con- 
struct physical densities from the densities l|6.10|l by way 
of 15.19(1 . The procedure is simply to replace e~ k / 4 with 
e -fc /2 -therein . We make the Fourier transformation of 
the physical densities with the aid of the formula [see 
formula (2.19.12.6) in Ref.j^] 



kdkk m e- k2,2 J m {kr)L { ^ (fc 2 /2) 



1 IT 



„2 \ n 



! V 2 



re 



-r 2 /2 



(6.12) 



and 



dO 



ikx 



2^J (fcr), 



J_ IdB^e ** = ^-Ji(fcr) 
2ni J k r n ' 



(6.13) 



in the coordinate space. 

In order to make a further analysis we make an 
anzats|l7j on functions u(m) and v(m), 



u 2 (n) = 



n + l + LO 2 



v 2 (n) = 



71+1 



n + 1 + u 1 



(6.15) 



The densities can be expressed in terms of the Kummcr 
function M(a; b; x), 



M(a;a+ l;i) = (iV 



n=0 



(n + a)n\ 



(6.16) 



P(x) 
S,(x) 

5j,(x) 



1 



1 1 



2?r 2tt + 1 



M(w 2 ;^ 2 + 2;r 2 /2), 



1 



1 



4ir 47T 
1 uj 2 



e-r ■ M(u 2 ;u 2 + l;r 2 /2) 



4tt w 2 + 1 



e - a p M(w 2 + l;w 2 + 2;r 2 /2), 



' ^ X -e-^ 2 A^ 2 + l; W 2 + 2;r 2 /2), 



Attlo 2 + 1 



^^ e ^M(c, 2 + l;. 2 + 2;r 2 /2), 

47T U) z + 1 



(6.17) 



where we have used the formula (13.4.3) of Rcf.^2J| to 
derive p(x). The electron number of the skyrmion exci- 
tation is 



= I d 2 x 



1 



e- z M(uj 2 ;uj 2 + 2;z)dz = -1, (6.18) 



where the last equality is checked via order by order inte- 
gration with respect to z. There is no ambiguity in this 
derivation contrary to that in l|6.11|) . 

The skyrmion with the anzats (|6.15() has a peculiar 
feature. It is reduced to a hole for uj — 0, where the 
density p(x) approaches the ground-state value exponen- 
tially fast. However, for all u> ^ 0, with the use of the 
formula (13.1.4) of Ref.jH], we find 



If 2a 2 \ 



as r — > oo, 



(6.19) 
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where we have set a 2 = 2w 2 . Furthermore, we find 



1 ax 



-1 ay 



Sz( x ) ~> 7, J> ^ O 2"' 



1/1 a" 



2tt V 2 r 2 



as r — > oo. 



(6.20) 



They agree with the asymptotic behaviors of the density 
and the spin field of a sufficiently large skyrmion we dis- 
cuss later: See IjS.lfill . It is easy to see that the number 
of the flipped spin diverges for all lo ^ 0. Consequently, 
small skyrmions cannot be discussed based on this an- 
zats. 



B. SU(4) Skyrmions 

The generalization to the SU(4) bilayer system is 
straightforward. The ground state is given by the up- 
spin bonding state. The skyrmion state is given by 



i6 Bky )=nf t (»)i°). 



(6.21) 



where 



£t( n ) = u ,(n)A\{n) + u p {n)A\(n) 

+ u s (n)B\(n) + v{n)B\{n + 1) (6.22) 

with constraint, u 2 .(n) + u 2 (n) + u 2 (n) + v 2 (n) = 1. The 
hole state £?f(0)|Bt) is given by setting v(n) = 1 and 
u s (n) = Up(n) = u T (n) = for all n. 
There are three types of antiskyrmions, 



|e s asky ) = nciW t (o)|o), 

n=0 

|SJU y >=II4(n)^t(0)|0) ) 

n=0 

i©L ky )=n^( n ) A i( o )i°)' ( 6 - 23 ) 



n=0 



with 



Ct(n) =u,(n)A\(n) +u p (n)A\(n) 

- w 8 (n)Bj(n + 1) + v{n)B\(n + 1), 

Cl(n) =Ur(n)A\(n) - u p (n + l)A\(n) 

+ u 8 (n)sj(n) +v(ri)B\(n + l), 

d(n) = - u r (n)A\(n + 1) + u p (n)A\(n) 

+ u s (n)B,(n) + v(n )B\(n + l). (6.24) 

They are reduced to three different electron excited states 
Sj(0)|BT), A|(0)|BT) and Aj(0)|Bt> when we set v(n) = 
1 and u s (n) — u p (n) = u T {n) = for all n: See FIG^ 

There are important SU(2) limits of SU(4) skyrmions. 
When we set u p (n) = u T (n) = 0, |e^ ky ) and |6| sky ) 



describe a skyrmion and an antiskyrmion where only 
spins are excited. We call them the spin-skyrmion and 
the spin-antiskyrmion. Similarly, when we set u s (n) — 
u r (n) = 0, |Sg ky ) and |©a sky ) describe a skyrmion and 
an antiskyrmion where only pseudospins are excited. We 
call them the ppin-skyrmion and the ppin-antiskyrmion. 

It is a dynamical problem which skyrmion-antiskymion 
pairs are excited thermally. As we have shown, as far as 
electron-hole pairs are concerned, only pseudospins are 
excited at the balanced point unless the tunneling gap is 
too large, while only spins are excited at the monolayer 
point. This is the case also for skyrmion-antiskyrmion 
pair excitations. However, in general, all components are 
excited to lower the total energy, which leads to genuine 
SU(4) skyrmions. Contrary to the case of the electron- 
hole limit, the transition from the ppin-skyrmion at the 
balanced point (ctq — 0) to the spin-skyrmion at the 
monolayer point (<7o = 1) will occur continuously via a 
genuine SU(4) skyrmion since the matrix elements of the 
total Hamiltonian H between various skyrmion states are 
nonvanishing. 



VII. EFFECTIVE HAMILTONIANS 

It is very hard to calculate the skyrmion excitation en- 
ergy with use of the microscopic states (|6.21f) and (|6.23|) 
since they involve infinitely many functions u s (n), u p (n), 
u r (n) and v(n). In this paper we construct the effective 
theory by making the derivative expansion of the Hamil- 
tonian. Thus, strictly speaking, our approximation is 
good only for large skyrmions. Nevertheless, its applica- 
tion even to small skyrmions would present us invaluable 
results otherwise unavailable. We wish to develop a mi- 
croscopic theory in a future work. 

Our analysis is based on the decomposition formula 
(|3.12|) into the direct and exchange energy terms. In what 
follows we represent the classical density p cl (x) simply by 
p(x) since we use only classical fields. 

We first identify the SU(4)-invariant direct Coulomb 
term as the self-energy, 



d 2 xd 2 yp(x)V + (x-y)p(y). 



(7.1) 



We then make the derivative expansion of the SU(4)- 
invariant exchange term H^. We rewrite it as 

H X=-\j d2xd " z ^x + (z) [I(x + z)T(x) 

+ ip(x + z)p(x)]. (7.2) 



Since Vx(z) is short ranged, 

^W=^ Jd 2 pe^Vi(p) 
=V(x) (l±e-^ d/e 



2 12l\ 



(7.3) 
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it is a good approximation to make the Taylor expansion 
of I(x + z) and /3(x + z) to the nontrivial lowest order of 

I(x + z)I(x) = T(x)T(x) -W z J z J 5 i I(x)5,i(x), (7.4) 



where a partial integration is understood in the integrand 
of 1)7.2)) . Equivalently, we make the momentum expansion 
of the potential 1)3.6)) . 

^x(p) = ^(0) - 8< J±p 2 + 0(p 4 ), (7.5) 

where 

K±(0) = 24 



1 ± e rf2 /24 erfc ^/V2£ B )1 Ac, (7.6) 



and 



Jf = ^{j s ±jf) 



with 



(7.7) 



(7.8) 



(7.9) 

The zeroth order term in p 2 is proportional to the 
integral 



of two planes with separation d, and the other (oc e x ) 
from the exchange interaction H^. The exchange ef- 
fect makes the capacitance energy quite small for a small 
layer separation 25j. We note that our capacitance for- 
mula 1)4.811 is different from the one assumed in some 
literature[6ll26l. 

Collecting all the first order terms in p 2 from the ex- 
change Hamiltonian 1)3.10(1 we obtain 



d k S a 



d k S a 



d k P a 



d k P z 



d k R a i 



d k Raz 



(7.13) 



where the summation over repeated indices a and b is 
understood. 

It is worthwhile to take two important limits of 7i x 
When all electrons are moved to the front layer, by set- 
ting W = (ip f \ V> , 0, 0), the nonvanishing elements are 
S a = Si and R az = S { a in f7~T3)> . and we find 



^spm 



2J, 

p\ 



d k Sl 



(7.14) 



for the spin-ferromagnet. Similarly, when the spin degree 
of freedom is frozen, the nonvanishing elements are P a — 
Pi and R az = Pj in 1(713) . and we find 



I(x)I(x) + ~?(x)p(x) 



1 



4< 



d 2 x /5(x), 



(7.10) 

which is a Casimir invariant obtained by integrating the 
relation 1)2.27)1 . Note that the star product becomes an 
ordinary product within the integrand. We may neglect 
the term, because it represents the total number of elec- 
trons and is fixed to the ground-state value in excitations 
of skyrmion-antiskyrmion pairs. Consequently, we obtain 



d fe T(x)d fc T(x) + —d k p(x)d k p(x) 



as the effective Hamiltonian. Here and hereafter the sum- 
mation over the repeated index in d k is understood. 

We next derive the effective Hamiltonian from the 
SU(4)-noninvariant terms Hq. The zeroth order term 
in p 2 yields the capacitance energy as the leading term, 



#cap = 27r4e cap / d 2 xP z (x)P,(x), 



(7.12) 



where we have used the identity 1)7.10)1 . The capaci- 
tance energy (oc e cap ) consists of two terms; the one 
(oc fp) arising from the direct interaction H^, which 
is the standard capacitance energy of a condenser made 



(7.15) 



for the pseudospin-ferromagnet. 

To discuss the Goldstone modes we may set p(x) = p 
since the ground state is robust against the density fluc- 
tuation. By setting S f a = p$,S a , 1)7.14)) becomes an 0(3) 
nonlinear sigma model describing the spin-ferromagnet 
Hamiltonian 7] with the spin stiffness J s . By setting 
Pa = Pi'Pa, 1)7.15)1 becomes an anisotropic 0(3) nonlin- 
ear sigma model describing the pseudospin-ferromagnet 
Hamiltonian[6j with the interlayer stiffness Jf. 



VIII. SEMICLASSICAL ANALYSIS 

We use bosonic variables to describe coherent excita- 
tions such as spin and pseudospin textures. In so do- 
ing we introduce the composite-boson (CB) field pll^. 
The CB theory of QH ferromagnets is formulated as 
follows^2- The CB field <^ CT (x) is defined by making a 
singular phase transformation to the electron field ip a (x) , 



ff (x) 



-ie6(x) 



^<t( x ). 



(8.1) 
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where the phase field 0(x) attaches one flux quantum to 
each electron via the relation, 

EijdidjQix) = $ D p(x). (8.2) 

We then introduce the normalized CB field n CT (x) by 

f(x)^(x)/(x), (8.3) 

where t (x)0(x) = E^lW^W = /°( x ), and 
the 4-component field 7i CT (x) obeys the constraint 
£ CT ra CT t( x )n ff (x) = 1: Such a field is the CP 3 field^. 
The isospin operators are expressed as 

5 (x) =p(x)5 (x), 5 a (x) = int(x)rf n n(x), 
P a (x) =p(x)P (x), P a (x) = int(x)rPP in n(x), (8.4) 

and so on, with n = (n 1 ^ , n 1 ^ , n h ^ , n b ^) T . They are 

1 



5| = 0, 5^ = 0, 5f= 2 



1 



1 



^I=2V 1_CT 0' ^ = ' V * = 2 a °> (8 ' 5) 

and ft* 6 = except for Tl\ x = V% and TZf z = Vf in the 
ground state. 

We investigate charged excitations at v — 1. Charged 
excitations are topological solitons in incompressible QH 
liquids. To describe them we introduce the dressed CB 
field by[l[l2 

^ Q (x) = e-" 4 ( x V Q (x), (8.6) 
where Q (x) is given by i|8.1[l and -A(x) is given by 



p(y) - M 2 - (8.7) 



A(x) = / dry In 



21 } 



The kinetic Hamiltonian (12.11) is rewritten as 



where $(x) is the four-component dressed CB field with 
$t = $t e 2^ 5 and £>. _ _ ifid . + h ( £jk + 5 jk )d k A(x). 

The LLL condition follows from the kinetic Hamilto- 
nian H8.8j) . 

(V x + l V y ) d>(x)|6) = A$ (x) | 6) = 0, (8.9) 

where z = (x + iy) /2is- It implies that the TV-body wave 
function is analytic and symmetric in N variables, 

6 C bH = (0|$(xi) • • • $( Xjv )|e). (8.10) 

It is easy to verifypj that the electron wave function 
is S[x] = 6cB[z]6LaughlinM, where 6 L au g hiin[x] is the 
Laughlin wave function j3JJ]. 



The analysis is quite simple when the function is fac- 
torized, 6cb[z] = Yl r &(z r ). Then it follows that &(z) — 
($(x)) and that (n(x)) = 6(z)/|6(z)|. The lightest 
topological soliton is described by the nontrivial simplest 
wave function (6 BT , & Bl , 6 AT , 6 Al ) = {z 
which we call the SU(4) skyrmion. Thus, a skyrmion 
is characterized by its shape parameters k s , k p and n r 
representing how it is excited to energy levels B|, A| 
and AJ,, respectively [FIG^. In terms of the layer CP 3 
field it reads 



/ n f T(x) \ 
n f i(x) 
n b T(x) 



= C(z) 



I Zy/ l + Q-Q + Kp\/ 1 — CTq ^ 
Zy/1 — OQ — K V \J\ + (To 



y K s \/l - (T - KVa/1 + O"0 / 



.11) 



with the normalization factor C(z) = 1/ \j2(z 2 + n 2 ) 
with k 2 — k 2 + k 2 + n 2 . When k = n s , k p — n r = 0, it is 
reduced to 



/n f T(x)\ 

n f i(x) 

n b T(x) 
\?i bi (x) I . 

V v ' / spi 



C(z) 



^ z yi + (TQ ^ 

Kt /1 + 

2: a /1 ^ (T^ 

\«vi - (t y 



.12) 



which describes a spin texture reversing only spins: This 
is identified with the microscopic spin-skyrmion in H6.21[l . 
When k = K p , Ks = n r = 0, it is reduced to 



/ n f T(x) \ 
n f *(x) 
n b T(x) 

V« bi (x)/ 



C(z) 



' Z a/1 + CT + Kx/1 - (T 





2 a/1 — 00 — Ka/1 + (To 

. 

(8.13) 

which describes a pseudospin texture reversing only pseu- 
dospins: This is identified with the microscopic ppin- 
skyrmion in l|tj.21f> . 

A skyrmion excitation modulates the density around 
it, p — > p sky (x), according to the soliton eauation[T2l|. 



— V 2 lnp sky (x) - p sky (x) +p = Jsky(x), 



;8.14) 

which follows from the LLL condition l|8.9|l : J s k y (x) is the 
topological (Pontryagin number) density, which is calcu- 
lated as 

4(^ B ) 2 



1 



Jsky(x) = - 



(8.15) 



7T [r 2 + 4( K e B ) 2 } 2 

for the SU(4) skyrmion (|8.11|l . The soliton equation is 
solved iteratively, and the first order term is 



^P s ky(x) ^ -^sky(x) 



1 



4( K £i 



(8.16) 



7T [r 2 +A(kI b ) 2 ] 2 ' 

This is good for a very smooth skyrmion (k> 1). 

The antiskyrmion configuration is related to the 
skyrmion configuration by 

^asky ^sky ^asky ^sky ^asky ^sky 

-pasky _ -psky -pasky _ _psky pasky _ psky 
' x 'xi ' y y 1 z ' z 

(8.17) 
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A skyrmion (antiskyrmion) induces the modulation of 
the spin and the pseudospin, 

*S«**(x) = p (a)sky (x)5( a ) sk "(x) - p S%, 
SP (^y {x) = P(8>ky (x)pW*(x) - p V*, (8.18) 

where and are the ground-state values (|8.5|l . and 



PskyW =P0 + < 5 Psky( X ); 
Pasky( x ) =Po - ^skyM, 



(8.19) 



where <5p sky (x) is given by l|8.16|l . 

It is actually the energy of a skyrmion-antiskyrmion 
pair, 



^pair 2 (^sky -^asky)> 



.20) 



that is observed experimentally. We estimate the exci- 
tation energy -E sky = (&\H\6) of one skvrmion[3l]. For 
simplicity we set n r = in the SU(4) skyrmion l|8.11|l . 
This approximation reduces the validity of some of our 
results. Indeed, we have found that the mixing between 
the excitation modes to the s-level and the r-level lowers 
the energy of the electron- hole pair state [FIG [2]- Never- 
theless, we use this approximation to reveal an essential 
physics of SU(4) skyrmions, because otherwise various 
formulas become too complicated to handle with. By 
an essential physics we mean a continuous transforma- 
tion of the SU(4) skyrmion from the ppin-skyrmion limit 
to the spin-skyrmion limit in contrast to the case of the 
electron-hole excitation. 

Thus we study a skyrmion parametrized by two shape 
parameters k s and k p with We calculate 

the SU(4) generators l|8.4|) for this field configuration, 



r 2 + a 2 



I a 2 



S sky = - 

2 r z + a 1 



ky = cos/? _ a p xsm(3 + a 2 p cos/3 = g p y 

2 r 2 +a 2 ' y r 2 +a 



sky _ 



sin P a p x cos (3 — a 2 sin (3 



.21) 



and similar expressions for 7?.^ y , where a s — 2k s £b, a p — 

2k p £b, a 2 = 4(k£b) 2 , cos/3 = \J\ — <j\ and sin/3 = <jq. 

The skyrmion energy consists of the Coulomb energy, 
the Zeeman energy and the pseudo-Zeeman energy. The 
Coulomb energy consists of the self-energy, the capaci- 
tance energy and the exchange energy. 

The dominant one is the self-energy l|7.1|l . 

£U = \f d 2 xd 2 y5 Psky (x)V+(x - y)Sp sky (y). (8.22) 
After a straightforward calculation we find 

E scli = j-E% J z 2 \K x {z)f (l + e-^) dz. (8.23) 



See (|C5|) in Appendix |0 It depends only on the total 



skyrmion scale k = y k 2 + k 2 , since the Coulomb term 

-ffp depends only on the total density 5p sky (x). 
The capacitance energy is given by (|7.12|) . or 

ff cap = 27r4e cap J d 2 x <5Pf y (x)<5P z sky (x). (8.24) 

It is calculated in (|C7fl in Appendix lO The leading term 
is 

-Bcap — ^(1 - crg)e C apA r ppi n (Kp), (8.25) 

where N pp i n (K p ) is the number of flipped pseudospins to 

be defined by fHH) . 

In evaluating the exchange energy (|7.13() . we set p(x) = 
p since the skyrmion charge is spread over a large do- 
main. (See Appendix B for the result without making 
the approximation.) Using (|8.21|) we obtain 



E x = 4tt 



' (8.26) 

It contains shape parameters explicitly via the SU(4)- 
noninvariant term. The SU(4)-invariant part of l|7.13|) is 
the SU(4) nonlinear sigma model yielding a topological 
invariant value AitJf . In the spin-skyrmion limit it is 
reduced to 



E s r = ^{jt+j^i) 



.27) 



which yields the well-known formulaQ, E^ m = 47rJ s , 
at the monolayer point (ao — * 1). In the ppin-skyrmion 
limit it is reduced to 

E^ = MJt - \jr) = f {Js + 2J s d ) , (8.28) 

which is independent of the imbalance parameter cto- 
The Zeeman energy of one skyrmion is 

E z = -A z Jd 2 x6S s J iy (x) = N spin (n s )A z , (8.29) 

where iV S pi n is the number of flipped spins. Neglecting the 
term which is cancelled out in a skyrmion-antiskyrmion 
pair excitation due to the relations (|8.18() . we obtain 



±{k s Ib) 2 
2 +4( K £ B y 



k 2 N{: (8.30) 



with 



Ns = 21n[l + -^-), (8.31) 

where the divergence has been cut off at r ~ k£ with a 
typical coherence length £. 
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The pseudo-Zeeman energy (|4.15|l is 



E 



PZ 



-A 



SAS 



d l x 



sp^y (x) 



<x &Ffr(x) 



(8.32) 



We extract the terms which are not cancelled out in a 
skyrmion-antiskyrmion pair excitation. It is equal to 



E 



PZ 



^ppin(^p) 



A 



SAS, 



where -/V pp i n is the number of flipped pseudospins, 

Nppin(Kp) ~ 

We have cut off the divergence as in l|8.3U[) . 
The total excitation energy is 



.33) 



.34) 



E s ky — E^ 



E, 



self - 



cap 



E z + E 



PZ 



(8.35) 



with f£U, fiTZty . and l|g35J) . Skyrmic 



parameters k s and re p with k 



k 2 , are to be 



determined by minimizing the excitation energy l|8.35[l . 
According to an examination of i|8.35|) presented in Ap- 
pendix^ provided 



4^ /1 

3 > U £ca P 



A, 



iSAS 



(8.36) 



ppin-skyrmions are excited at the balanced point (<ro = 
0): See <|D3|) . Then, it evolves continuously into a spin- 
skyrmion at the monolayer point (co = 1) via a generic 
skyrmion {k s k p ^ 0). As we have stated, to simplify 
calculations we have set n r = in the SU(4) skyrmion 
(|8.11|l . When we allow a skyrmion to be excited into the 
r-level (n r ^ 0), a genuine SU(4) skyrmions (K s K p K r ^ 
0) would be excited except at the monolayer point, as 
illustrated in FIGdl Compare this with FIG0 

We examine the condition i|8.36[l numerically. When 
we adopt typical values of sample parameters (p — 
1.2 x 10 11 /cm 2 and d = 231 A), we find the capacitance 
(|4.8|) to be e cap ~ 132K while the exchange-energy dif- 
ference to be 47rJ~/3 ~ 5.4K. The condition is hardly 
satisfied. Here, we question the validity of the standard 
identification of the layer separation, d = ds +dw, where 
dw — 200A is the width of a quantum well and ds — 3lA 
is the separation of the two quantum wells. In this identi- 
fication the electron cloud is assumed be localized in the 
center of each quantum well. However, it is a dynamical 
problem. Let us minimize the ppin-skyrmion energy as a 
function of d. The minimum is found to be achieved at 
d < ds- Namely, the energy increases monotonously for 
d > ds- Then, it would be reasonable to use d = ds as 
the layer separation to estimate the energy of the ppin- 
skyrmion. When we choose d ~ 3lA we find e cap — 4.5K. 
Then, the condition i|8.36[l is satisfied in some parameter 
regions, where ppin-skyrmions are excited. 

However, the condition i|8.36|) cannot be taken literally, 
because it rules out excitations of large ppin-skyrmions 
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FIG. 3: The skyrmion energy is estimated as a function of the 
imbalanced parameter <to with sample parameters Asas = 
IK, and p — 1.2 x 10 11 /cm 2 . The experimental data are 
taken from Terasawa et al.jTg]. Solid lines show the energies of 
a ppin-skyrmion and a spin-skyrmion, normalized to the data 
at the balanced point and the monolayer point, respectively. 
The experimental data would be explained by the excitation 
of a genuine SU(4) skyrmion. 



(k 2 N^ 3> 1) against experimental indications [Til fbll IT U . 
It seems that the exchange energy has been underesti- 
mated by making the derivative expansion. It is nec- 
essary to go beyond the present approximation to fully 
understand the problem. 

We proceed to estimate the energies of a spin-skyrmion 
and a ppin-skyrmion in all range of tro, because their es- 
timation is much more reliable than that of a generic 
SU(4) skyrmion. We use sample parameters p = 1.2 x 
10 11 /cm 2 , d = 3lA and Asas = IK. and normalize the 
energy to the experimental data(3] a t the monolayer 
point for a spin-skyrmion and at the balanced point for 
a ppin-skyrmion. Recall that the absolute value of the 
activation energy cannot be determined theoretically be- 
cause they depends essentially on samples. 

First, we calculate the energy of a spin-skyrmion as a 
function of (To by setting k = k s and k p — in l|8.35|l . 
The skyrmion scale k, to be determined by minimizing 
this, depends on ao very weakly. We obtain N sp i n (K) ~ 
3.5 as in the monolayer QH ferromagnet [j, [l2j . We have 
depicted the excitation energy in FIGEH where it 

is normalized to the data at the monolayer point. 

Next, we calculate the energy of a ppin-skyrmion as a 
function of ao by setting k = k p and k s — in l|8.35|l . 
There is an important remark. The effective tunneling 
gap diverges as Co — > 1 unless Asas = 0. Further- 

= 1 for a 
Hence, the pair excitation 



more, the minimum pseudospin flip is iV, 
skyrmion-antiskyrmion pair, 
energy l|8.20|l diverges as 



iSAS 



— ■ 



(8.37) 



and ppin-skyrmions are not excited at the monolayer 
point (ctq ~~ * !)• We have depicted the excitation energy 
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-^sky 11 m FIG0 where it is normalized to the data at the 
monolayer point. In so doing we have used the full ex- 
pression l|C7(l for the capacitance energy E cap since other 
terms are as important as <|8.25[1 for a ppin-skyrmion of 
an ordinary size. 



IX. ACTIVATION ENERGY ANOMALY 

We have studied how one skyrmion evolves continu- 
ously from the balanced point (co = 0) to the monolayer 
point (cto = 1) by changing its shape. It is important 
how to distinguish various shapes of skyrmions experi- 
mentally. As is well known, as the sample is tilted, the 
activation energy of a spin-skyrmion increases due to the 
Zeeman energy. On the contrary, the activation energy 
of a ppin-skyrmion decreases due to the loss of the ex- 
change energy. Thus, the tilted-field method provides us 
with a remarkable experimental method|lllll4| to reveal 
the existence of various shapes of a skyrmion in bilayer 
QH systems. 

As the sample is tilted, the parallel magnetic field Bn 
is penetrated between the two layers. We take the sym- 
metric gauge generalized as 



A = 



(9.1) 



where the two layers are placed at z = ±g?/2. In the ki- 
netic Hamiltonian (|2.1|) the covariant momentum is now 
different between the two layers, 



D T = — ih 



■ t d 


h 


th 7T + 
ox 


W B 


■n 9 


h 


dy 





l) 



eB\\d 



(9.2) 



Consequently, the LLL condition (|8.9(l is modified as 



{V x +iV y )^)\e-B^ 
ih 



with 



|;+^jBrr) ( I > (xi!6:// i 



edBu 





(9.3) 
(9.4) 



We may solve the LLL condition (|9.3() for the one-body 
wave function 6(x;_Bii) = (x|6;£?ii) as 



&( x ; B \\) = ex P 



6 m T p z pm x I 6(z). (9.5) 



Accordingly, the skyrmion configuration acquires differ- 
ent phase factors between the two layers, 



n(x;i?||) = cxp( --<5 m rPP in ^)n(x;0), (9.6) 
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FIG. 4: The activation energy at v = 1 is plotted as a func- 
tion of the tilting angle in several samples with different 
tunneling gaps Asas- It shows a rapid decrease towards the 
critical angle O*, and then becomes almost flat. The data 
are taken from Murphy et al.|lll|. They are well fitted by 
the theoretical formula 19. lit , where the activation energy 
decreases in the commensurate phase and becomes flat in the 
incommensurate phase. To fit the data, we have adjusted 
the activation energy at O = with the experimental value, 
and assumed that the number of flipped pseudospins N vvm is 
constant for all values of the tilting angle. 



where n(x; 0) is the configuration (|8.11|l in the absence 
of the parallel magnetic field. Various isospin fields are 
given by with this CP 3 field. 

We first consider the balanced configuration (er = 0), 
where we assume ppin-skyrmions are excited [FIG OH- 
The excitation energy at B« = is given by, 



£ s p k p ; n (o) 



£P pm (0) 



E f 



olf ■ 



Nppi n A S AS , 



(9.7) 



where E§? ln 



(0) is the exchange energy (I8.28|) in the ab- 
sence of the parallel magnetic field. We analyze how the 
exchange Hamiltonian (|7.15|l is affected by the parallel 
magnetic field. Subtracting the ground-state energy we 
easily deduce the Bu dependence of the excitation energy, 



££ pin (B||) = £P pin (0) + A£P pin (B||) (9.8) 



with 



AE^(B l{ ) = -2jf6 2 m fx (^(x)) 



(9.9) 



where "Pf ky is the skyrmion pseudospin component in the 
absence of the parallel magnetic field and given by H8.21J) . 

We note that AE^ ln (B\\) is proportional to the ca- 
pacitance energy l|7.12|l . Thus the leading order term is 
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found to be 



2TTd 2 J° s 

I 2 

B 



iVppi n tan 2 6, 



We proceed to analyze excitations of SU(4) skyrmions 
(|9.6|) with H8.11fl iii imbalanced configuration, where 
(9 10) the exchange Hamiltonian is given by (|7. 1 . After 
some calculation we find that the exchange-energy loss 
is again proportional to the capacitance en- 



where tan 6 = B»/B± and A DDin is the number of pseu- x , , , J 1 , ,. , 

, , , , . , ergy, and the leading order term is 



dospins flipped around the skyrmion. All other terms 
in (|9.7(l are unaffected by the parallel magnetic fieldQ. 
Thus, the excitation energy decreases as the tilting angle 
© increases. The rate of the decrease depends on the 
number A PP i n of flipped pseudospins and the amount of 
the penetrated magnetic field Bu . 

It has been shown[l|, |6j that, when the parallel mag- 
netic field increases more than a certain critical point, 
the phase transition occurs in the bilayer QH system: It is 
the commensurate- incommensurate transition point B^ . 
In the incommensurate phase the penetrated magnetic 
field is not increased more than BT,, because the excess 
magnetic field A-Bjj = £?* — BT, is eaten up to create pen- 
etrated sine-Gordon vortices between the layers. Hence, 
the activation energy becomes flat for > 0*, where 
tan©* = BpB ± . 

Hence, from (|9.8ll and l|9.10|l the excitation energy is 



-E, 



self" 



-E, 



cap 



-Ar ppin Af AS , (9.11) 



where we have absorbed AE^ m (B^) into the effective 
tunneling gap, 



i SAS 



Asas 

AsAS 



(2ird 2 jf/£ 
{2vd 2 J d Jl 



tan 2 for < 0* 

tan 2 0* for > 0* 

(9.12) 

At each angle 0, in principle, it is necessary to min- 
imize the total excitation energy H9.11(l and determine 
the scale parameter k. However, since the energy change 
(|9.12l) is quite small compared with the total energy, we 
may treat it as a small perturbation. Namely, we may 
assume that the flipped pseudospin number A pp i n is a 
constant independent of 0. This reminds us of experi- 
mental data 32, 33], where the flipped spin number N sp i n 
is nearly constant by tilting the sample in the excitation 
of spin-skyrmions. 

We have fitted the data due to Murphy et al.[ll| in 
FIG 01 by assuming an appropriate flipped pseudospin 
number A pp j n per one skyrmion-antiskyrmion pair. We 
summarize the result in a table. 



sample 


A 


B 


C 


D 


Asas 


0.81 


0.86 


4.43 


8.53 


A PP in 


18 


24 


2.1 


1.5 



(9.13) 



It takes a large value in samples with Asas < 1 K but 
take a small value in samples with Asas > 4 K. Recall 
spin-skyrmion excitations in the monolayer QH system, 
where the flipped spin number remains small when the 
Zeeman energy is moderate but becomes quite large when 
the Zeeman energy is almost zero. 



A£; x (i3||)^-^^(l-a 2 )A ppin tan 2 0, (9.14) 

B 

which is reduced to (|9.10|l in the balanced point. 

Hence, in the commensurate phase (0 < 0*) the ex- 
citation energy turns out to be 



£ sky (5||) =E X (0) + E sclf + JScap+JVgpiaVl + tan 2 9A Z 
+ A ppin Af AB (9.15) 

with 

1 . 2ird 2 jf 



^BAB 



^SAS~ 



f 2 

B 



-(l-o-^tai^e. (9.16) 



In the incommensurate phase (0 > 0*) it is given by 
this formula by replacing tan 2 © with tan 2 0*. Here, 
the commensurate-incommensurate transition point in- 
creases slowly as[34| 



(i-^r 1/4 e* 



(9.17) 



as the imbalance parameter ao increases. 

Experiments have been carried out by Sawada et al.|l4| 
in bilayer samples (FIGEJ), where the activation energy 
was measured at v — 1 by controlling both the tilting 
angle and the imbalance parameter cto- Their data 
are interpreted based on the theoretical result (|9.15(l as 
follows. We focus on the behavior of the activation energy 
-E^ky by changing the tilting angle at each fixed imbalance 
parameter. The relevant term in l|9.15ll is 



A£ skv (0) =A spin v/l + tan 2 OA° z 

2^(1 - a 2 ) 



ppm 



( 2 
' B 



tan 2 © (9.18) 



for © < ©*. By adjusting the theoretical curve to the 
data at the point = 0, we fit the data by this curve 
throughout the observed range of the tilting angle 0. As 
seen in FIG|5l the fitting is quite good when by assuming 
constant values of A sp i n and A pp j n throughout the range 
of 0. A deviation of the theoretical curve from the data 
for large tilting angles > 70° would be due to effects 
not taken into account in the above analysis. For in- 
stance, when the parallel magnetic field become too large, 
the soliton lattice becomes too dense in the incommen- 
surate phase and would destabilize skyrmions. We sum- 
marize the numbers A sp ; n and A pp i n per one skyrmion- 
antiskyrmion pair determined by this fitting in a table. 
As the sample is tilted, A sp j n increases and A pp j n de- 
creases. In particular, both spins and pseudospins are 
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FIG. 5: The activation energy at v = 1 is plotted as a func- 
tion of G in one sample (Asas = 11 K, p = 0.6 x 10 11 /cm 2 ) 
with different imbalance parameter ao. It shows a decrease 
towards the critical angle 0*, and then begin to increase for 
o"o = 0.3 and ao = 0.7. The data are taken from Sawada et 
al.|l4j|. To fit the data by the theoretical formula (I9.15H . we 
have adjusted the activation energy at O = with the ex- 
perimental value, and assumed that the nipped spin number 
iVppi n and the flipped pseudospin number N pp i n are constant 
for all values of the tilting angle. 

flipped unless <7n = or (To = 1. We conclude that 
the SU(4) skyrmion evolves continuously from the ppin- 
skyrmion limit to the spin-skyrmion limit as un increases 
from to 1. 








0.3 


0.6 


0.7 


1 


N~ spin 





0.64 


2.3 


2.8 


7.0 


■*V; ppin 


1.9 


1.8 


1.5 


0.68 






It is interesting to study the same problem in a sample 
having a very large tunneling gap. Terasawa et al.|l5j 
have measured the activation energy by controlling the 
tilting angle O at the balanced point in a sample with 
A S as - 33K. We have fitted their data [FIG® by the 
ppin-excitation formula l|9.11f) and by the generic formula 
(|9.15|) . It is difficult to fit the data if pure pseudospin ex- 
citations are assumed since it is required that -/Vp P i n > 1 
per one pair. A better fitting is obtained if spins and 
pseudospins are excited simulaneously since N pp i n may 
take a smaller value than 1. Such a simultaneous exci- 
tation is allowed at the balanced point as explained in 
Appendix iDl 

In passing we comment on the original mechanism 6, 
[35|,|36j proposed to explain the activation energy anomaly 
based on the exchange-energy loss of bimeron excitations. 
A bimeron has the same quantum numbers as a skyrmion, 
and it can be viewed as a deformed skyrmion with two 




M 40 60 9b 

Tilting Angle (degree) 



FIG. 6: The activation energy at v = 1 is plotted as a func- 
tion of in a sample (Asas = 33K, p = 1.0 x 10 11 /cm 2 ) . 
The data are taken from Terasawa et al.|l5j|. To fit the data 
by the theoretical formula, we have adjusted the activation 
energy at O = with the experimental value. We have also 
assumed that iV S pj n and N pp i n are constant for all values of 
the tilting angle. A better fitting is obtained as indicated by 
a solid line when both spins and speudospins are excited. 

meron cores with a string between them. The bimeron 
excitation energy consists of the core energy, the string 
energy and the Coulomb repulsive energy between the 
two cores. It is argued that the parallel magnetic field 
decreases the string tension and hence the bimeron ac- 
tivation energy. Clearly the mechanism works well only 
when the string length is much larger than the core size. 
A microscopic calculation has already revealed [37| that 
the meron core size is large enough to invalidate the naive 
picture. Furthermore, the skyrmion is almost as small as 
the hole itself in samples with large tunneling gap. On 
the contrary, in our mechanism the decrease of the ex- 
change energy follows simply from the phase difference 
induced by the parallel magnetic field between the wave 
functions associated with the two layers, and it is valid 
even for small skyrmions. 



X. DISCUSSION 

We have investigated the dynamics of bilayer QH sys- 
tems based on an algebraic method inherent to the non- 
commutative plane with [X, Y] — —i(? B . The noncom- 
mutativity induced by the LLL projection implies that 
the electron position cannot be localized to a point but 
to a Landau site occupying an area 2ir£ B . We have de- 
rived the Landau-site Hamiltonian Hq akin to the lat- 
tice Hamiltonian. It has two entirely different forms, the 
direct-interaction form Hy> and the exchange-interaction 
form Hx- They are equivalent, ifo = Hx, as the micro- 
scopic Hamiltonian. Nevertheless, the energy of a charge 
excitation consists of two well-separated pieces, the direct 
energy H$ and the exchange energy H x l . 

One of our new contributions is the derivation of vari- 
ous LLL-projected Coulomb potentials in analytic forms. 
For instance, we have revealed a new form of the capac- 
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itance energy (|7.12|) with e cap = 4(6^ — e^). It is a sum 
of the direct and exchange Coulomb interaction effects, 
where the direct effect yields the well-known result from 
the planar condenser proportional to the layer separa- 
tion, while the exchange effect gives a quite large nega- 
tive contribution to it when the layer separation is small 
enough. 

We have explored various aspects of SU(4) skyrmions 
at v = 1. In particular we have studied a skyrmion- 
antiskyrmion pair in its small size limit (electron-hole 
pair) and its large size limit. 

The excitation energy of an electron-hole pair is exactly 
calculable. We have obtained the excitation energy as a 
function of the imbalance parameter oq. The result is 
quite interesting: At the balanced point (<jq = 0) the 
pseudospin excitation occurs provided the tunneling gap 
is not too large (Asas — Az < 2e^ ). A peculiar feature is 
that, as <7o increases, the spin excitation occurs suddenly 
because the excitation involves just one electron or hole. 
At the monolayer point (cto = 1) °uly spins are excited 
always. We have also extended the microscopic theory of 
skyrmions |17| to our framework. However, a quantitative 
analysis is yet to be carried out. 

We then have estimated the SU(4) skyrmion excitation 
energy as a function of cfq based on the effective Hamil- 
tonian valid for very smooth isospin textures. In typical 
samples it flips only pseudospins at the balanced point. 
As (To increases, it evolves continuously to flip both spins 
and pseudospins, and finally flips only spins at the mono- 
layer point. We have then calculated how the excitation 
energy changes as the sample is tilted. Our formula has 
explained quite well the activation energy anomaly found 
by Murphy et al.[llj at the balanced point by excitations 
of pseudospins, and also found by Sawada et al.0] at 
various values of ao by simultaneous excitations of spins 
and pseudospins. Though our formulas are derived for 
sufficiently smooth skyrmions, they have turned out to be 
quite good at least with the use of phenomenological val- 
ues of iV S pi n and -/Vp P i n . We wish to develop a microscopic 
theory to analyze small skyrmions in a future work. In 
conclusion, the activation energy anomaly is explained 
by the loss of the exchange energy of SU(4) skyrmions, 
which are reduced always to spin-skyrmions at the mono- 
layer point and mostly to pseudospin-skyrmions at the 
balanced point. 



APPENDIX A: GROUP SU(4) 

The special unitary group SU(N) has (A^ 2 — 1) genera- 
tors. In the standard representation!^, we denote them 
as A^, A = 1, 2, • • • , N 2 — 1, and normalize them as 



Tt:(\a\b) = 2$ab- 



(Al) 



They are characterized by 
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{A,4,Ab}=— + 8 AB 2d A BC^c, (A2) 



where Jabc an d cIabc are the structure constants of 
SU(N). We have \a = ta (the Pauli matrix) with 
J abc = £abc and cIabc = in the case of SU(2). 

This standard representation is not convenient for 
our purpose because the spin group is SU(2)®SU(2) 
in the bilayer electron system with the four-component 
electron field as * = (i/) fT , ip^, i/> bT , 4> hl ). Embedding 
SU(2)<g>SU(2) into SU(4) we define the spin matrix by 



t v 

T„ 



spin 




(A3) 



and similarly the pseudospin matrix by 



ppin 




ppm 




(A4) 



where I2 is the unit matrix in two dimensions. Nine re- 
maining matrices are products of the spin and pseudospin 
matrices: 
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related with the standard SU(4) generators as 



provided 





— Ai + A13, T y o — A2 + 


A14, 








T zQ 


1 y/6 
^A 3 __A 8 + _ A 15 , 












= A4 + An, Toy = A5 + 
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Tqz 


= 7! A8 + 7i Al5 ' 
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= A6 + Ag, T xy = A10 + 
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Txz 


= Ai 


" A13, 
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= Mo — -V j Tyy = Ag — 
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= A 2 
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T zz 


-A + 1 A ^A 
_A 3 + - 7 =A 8 -— Ai 5 . 








(A6) 



The normalization condition reads 

Tr{T ab T cd ) = 4S ac S bd , (A7) 
which is different from the standard one ijAlfl . 

APPENDIX B: DECOMPOSITION FORMULA 

The Landau-site Hamiltonian Hq possesses two en- 
tirely different forms, the direct-interaction form Hd and 
the exchange-interaction form Hx- They are equivalent, 
Hd = -ffx, as the microscopic Hamiltonian. On the other 
hand, the excitation energy consists of two well-separated 
pieces, the direct energy and the exchange energy. In this 
appendix we derive the decomposition formula l|3.12[l . 

(&\Hc\&)=H$ + H* (Bl) 

for skyrmion excitations. We also prove the algebraic 
relation $2.27} . or 

N 2 — 1 

£ /l 1 (x)*71 1 (x) + J-p cl (x)*p cl (x) = -^(x), 



A=l 



(B2) 

that holds among the classical densities associated with 
the generators of the Woo(N) algebra. 

We consider a skyrmion state in the SU(N) QH 
fcrromagnet, 



\e) = HzHn)\o), 



(B3) 



n=0 



where 



N 



= £ + Mn)4(n + 1)] . (B4) 



They satisfy the standard canonical commutation 
relations, 

K(m),C + (n)] = S mn , K(m),C(n)] = 0, (B5) 



N 



[u»(n)ul(n) + Vp(n)vl{n)] =1, 

N 

5>t(n)u p (n+l)=0. (B6) 

A" 

It follows that 

C/4 (n)|6) = Mn)£(n) + v^n - l)£(n - 1)] |6), (B7) 

from which the only nonvanishing components of two- 
point correlation functions are found to be 

(c*(n)c„(n)) = u*(n)u„(n) + u*(n - l)v u (n - 1), 
(c* fl (n)c„(n + l))=u;(n)v v (n), (B8) 

where (c*(n)c u (n)) = (&\c*(ri)c v (ri)\&) and so on. We 
also derive 

c^j)c„(n)\&) = MJ)£(J) + v^U - IW - 1)] 

x [u„(n)f(n) + tv(n - l)f(n - 1)] |6). 

(B9) 

Four-point correlation functions (c^ (m)ct (i)c T (j)c u (n)) 
are complicated. Nevertheless, taking into account the 
angular-momentum conservation (V mn ij oc <5 m +i.„+j) we 
deduce 

V mrl ij (c^{m)cl{i)c T {j)c v {n)) 

=V mni j (cl(m)c u (n)) (4(i)c r (j)) 

- Vmnij <4(m)c T (i)) (4(i)c„(n)) . (BIO) 

This is the general expression allowing to express all kind 
of Coulomb energies in the exact form. 
It follows from (iBlOl) that 



Vmnij Y ( C jt( m ) C 3r(*) C r0') C ''( n )) 6 MV S aT 



cl/ 



=V mnij p c (m,n)p c - V mnij X(m,n,i,j), (Bll) 
where we have set 
X(m,n,i,j) = Y ( c t( m ) c T(j)) (4(*) c 4 n )) 8pu>6<rr- 

(B12) 



\1VTV 



Here we use the algebraic identity (|3.1[) to find 

1 



X(m,n,i,j) = 2 
Consequently, 



I%(m,j)I%(i,n) 



2N 



P (m,j)p cl (i,n) 
(B13)' 



mnij <y -J 

=E U + H x , 



(B14) 



20 



where 

Eb = Y. V mm]P c \m,n)p cl (i,j), (B15) 

innij 

E X = - 2 J2 V mmj [ J2 7 ^ m < 3)lA{i, n) 

mnij A—l 

+ ^p cl (m,j)p c \i,n)}. (B16) 

Now, Ey> and Ex are identical to (|2.18|l and l|3.2[l when 
we replace p(m,n) and J^(m, n) with p cl (m,n) and 
Ijl(m,j). Hence, we have established the decomposition 
formula l)Bljl rigorously. 

We proceed to prove the algebraic relation (|B2(1 for the 
skyrmion state (|B3|) . From l|B12|l we define 



X(m, n) = X(m, n, i, j)5ij 

= E E ( c U m ) c °(j)) <4o>») • ( B17 ) 



Using (|B6|) and (|B8|I we find the only nonvanishing values 
to be given by 

X(n, n) =v) (n)u(ri) + v\n — l)v(n — 1), 
X(n, n + 1) =[X(n + 1, n)]* = u\n)v{n), (B18) 

which implies that X(m, n) — p cl (m, n). Combining this 
with IB 1311 we obtain 



I c 2(m,l)F2(hn) + ± p °\ m ,l)p°% n ) = \p c \m,n). 

(B19) 



This is equivalent to (|B2|) in the coordinate space. 



APPENDIX C: COULOMB ENERGIES 

We estimate the Coulomb self-energy of one skyrmion 
in the bilayer QH ferromagnet. It is defined by 

E sclf = ir J d 2 qSp sky (-q)V + (q)Sp aky (q) (CI) 



with 



(C2) 



The density modulation i5p sky (x) is given by l|8.16[l for a 
large skyrmion. Its Fourier component reads 

WO) = ^J rf 2 ^e iqx ^ skv (x) = ^K x {aq), (C3) 

where we have set a = 2kIb- To calculate the energy we 
use (|6.13|) and 



* (r 2 +Q 2 )m+ i = -^T K ^{aq), (C4) 



where K n (z) — K_ n (z) is the modified Bessel function: 
See the formula (11.4.44) in Ref.[24|. Now it is straight- 
forward to derive 

E seli = |f / AKi(z)} 2 (l + e-**) dz (C5) 



as the Coulomb self-energy of one skyrmion. 

The SU(4)-noninvariant term yields the capacitance 
energy l|8.24ll . 

E cap = 2^|e cap J d 2 x 6P* y (x)5Pf y (x), (C6) 

where <5P| ky (x) is given by l|8.18|l . After some calculation 
we obtain 



E a 2 



c cap 



+ (1 - ai)K* p A{K) + ai v ^ 2 . . ^ 2 



1 



1 



(07) 



where a = 2k^b and a p = 2k p £b- A(k) is divergent 
logarithmically, 



2a 2 f 



2*2 



+ a 2 (r 2 +a 2 y 



5 dr. (C8) 



We take the divergent term as the leading contribution 
for the pseudospin component, 



1 



(.'■ a i > — (1 cr o)^ ca P'^ /r ppi n (^p) 



(C9) 



with ipOijl . 

Finally we present the formula for the exchange en- 
ergy l|7.13|) for a general SU(4) skyrmion (|8.11|) . After a 
straightforward but tedious calculation we obtain 



SU(4) 



E 



=4ttJ+ 1 



10k 4 



- 4tt Jr - 



1 l - 3af 



3 140k 4 



,'2 3 
4rr J7 - 



3 14k 4 / V K V K / K 

,2 1 in„2\ „2 



\ 2 .2 



An J: 



+ 4mJ-{l-oi) - 



1 + o-g l-19crg\ k; 



3 ' 140k 4 / k 2 

I I \ K 2 



+ 4ttJ s o- \/l- o\ 



3 ' 140k 4 / k 2 

2 ^ ^^7" 



3 70k 4 y k : 



(CIO) 



By setting n r = 0, this is reduced to E X V ^ = Ex + AEx 
with l|8.2(j[) . where 

AEx "lo^ Js ~ Ho^ J * 

+ ii^[ 3k4 - 19k2 ^+ 30 ^ j - ( cu ) 

The correction term is small for a large skyrmion. 
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APPENDIX D: CONTINUOUS 
TRANSFORMATION 

Based on the skyrmion-energy formula (|8.35|) we ver- 
ify that, if a ppin-skyrmion (k s = 0, k p = k) is ex- 
cited at the balanced point, it evolves continuously into 
a spin-skyrmion (re s = k, k p = 0) via a generic skyrmion 
(k s k p 0) as the imbalance parameter cro increases. 

We summarize the total skyrmion energy 1)8. 35|) as a 

function of k = \J k 2 + k 2 , z — (k s /k) 2 and (To hi the 
following form, 

4 7r j+ -2 

E Bky (K,z;a ) = j—9-z 2 +A(n;<j )z+B(n;ao), (Dl) 



where 



AttJ; 
T 



A z k 2 A^ 



A. 



SAS 2 



(D2) 



with (|8.31|1 . The explicit expression of B(n:ao) is not 
necessary. The variable z is limited in the range < z < 
1. 

We start with the balanced point (cro = 0), where 
^sky(^, 2; Co) is a linear function of z. Let us assume 



j4(k;cto)>0 at do = 0. 



(D3) 



Then the energy is minimized at z = or k s = 0, where 
ppin-skyrmions are excited. We now increases ctq from 



(To = 0. As far as A(k; <tq) > 0, ppin-skyrmions have the 
lowest energy. Due to the Asas term, A(k\<jq) decreases 
and vanishes at ctq = °~p- 



A(k; Co) =0 at (Tq 



(D4) 



Then it becomes negative for <to > er p , and the energy is 
minimized at z = z mul with 



8ttJ. 



\A(K;a )\. 



(D5) 



Because z m - ul increases continuously from z m - m — at 
°"o = ^p; the spin component is excited gradually 
to form a generic skyrmion (k s k p ^= 0). The point z m - ul 
increases and achieve at z m j n = 1 at ctq = <J S , 



\A(k;<t )\ 



8nJ+a 2 



at (T = (T s 



(D6) 



where the pseudospin component (n p ) vanishes continu- 
ously. 

We conclude that, if A(n; (To) > at (To = 0, ppin- 
skyrmions are excited for < (To < ^pj generic skyrmions 
are excited for a p < er < cr s , and finally spin-skyrmions 
are excited for a s < cro < 1. The transition occurs contin- 
uously as illustrated in FIG|3 with the critical points a p 
and (T s being fixed by (|D3J) and l|D6|l . On the other hand, 
generic skyrmions are excited at the balanced point if 
A(k; (To) < and z m - ul < 1 at <tq = 0, and spin skyrmions 
are excited at the balanced point if A(k; ctq) < and 
Zmin > 1 at cr = 0. 
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